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ABSTRACT 

Equilibrium configurations of weakly interacting fully degenerate fermionic 
dark matter are considered at various scales in the Universe. We treat the gen- 
eral situations for the gravity from Newtonian to general relativity and the de- 
generacy from nonrelativistic to relativistic. The formulation of the problem is 
exactly the same as the case treated by Oppenheimer and Volkoff in their paper 
on neutron stellar cores. A dimensionless equilibrium configuration is specified 
by a single parameter regardless of particle properties, the Fermi velocity at the 
center, and the scalings of mass and length are specified by the rest mass and 
statistical weight of the dark matter particle. We focus our attention to the flat- 
top nature of the mass column density profile of the cluster of galaxies, A1689, 
recently reported by Broadhurst et al. using gravitational lensing. We convert 
the column density profile to a volume density profile assuming spherical sym- 
metry and derive a 3D encircled mass profile of A1689, which is compared with 
the model profiles of degenerate fermion structures. The flat-top profile is repro- 
duced. The corresponding fermion mass ranges from 2 eV to 30 eV depending 
on the actual scale of the degenerate structure. If massive neutrinos are the 
dominant dark matter, the rest mass will be about 4.7 or 2.3 eV respectively 
for Majorana or Dirac neutrinos. The mass and size of the degenerate structure 
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are 1O 14 M and 100 kpc for Majorana neutrinos, and 5xl0 14 M Q and 300 kpc 
for Dirac neutrinos. If we identify the fermions as heavier sterile neutrinos, they 
yield the characteristic mass hierarchy of black holes; giant black hole at the 
center of a galaxy and the intermediate mass black holes. Thus we propose the 
possibility that the mass hierarchy of fermions determines that of black holes in 
the Universe. 

Subject headings: clusters:individual(A1689) — gravitational lensing — neutrinos 
— black hole physics 

1. Introduction 

Recent precision observations have revealed that the unknown dark matter dominates 
the matter contents of the Universe. We wish to study the possible dynamical structures of 
their existence, especially in a universal form. If the dark matter is in the form of ordinary 
thermal gas, the structure and the dimension would be strongly dependent on the initial 
conditions and the environment of the expanding universe. On the other hand if the dark 
matter is almost degenerate, we will naturally find a universal structure. Moreover, we would 
like to know the possibility that the dark matter forms black holes. 

To make the problem setting better defined, we start our consideration from the typical 
structure made from ordinary baryonic matter, in which an electron and a nucleon form 
the basic ingredient through the electromagnetic force. In the total energy E = \m e v 2 — 
(e 2 /r), the second term represents the attractive force making the system collapse and the 
first term represents the pressure against it through the Heisenberg's uncertainty principle. 
Actually putting the expression of the de Broglie wave length r = h/ (m e v) in E, and 
extremizing it with respect to r, we obtain the ground state bounding energy of hydrogen 
E = -m e e A /(2h 2 ) = -2.19 x 10~ 18 J and the Bohr radius r = H 2 /(e 2 m e ) = 5.28 x 10~ 9 
cm. The corresponding energy density is p = m p /(^-r 3 ) = 2.71 g-cm~ 3 . This is the basic 
ingredient of the structure formed by electromagnetism. The electron can collapse more. 
Setting v — > c in de Broglie wave length yields the Compton wave length. The corresponding 
energy density becomes p = rn p /(^-r 3 ) = 6.93 x 10 6 g-cm" 3 . This is the basic ingredient of 
white dwarfs. Further collapse makes the electron and the proton into a neutron, through 
the inverse beta decay process. Putting m e — > m p , we have r = h/(m p c), and therefore 
the corresponding energy density is p = m A v (?jY? = 4.29 x 10 16 g-cm~ 3 . This is the basic 
ingredient of neutron stars. 

In general, the material formed by fully degenerate fermion of mass m would have the 
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energy density p = m 4 c 3 / ft 3 . A structure of this density with radius R has the mass M = R 3 p. 
On the other hand the size R should be smaller than the limiting scale, Schwarzschild radius, 
2GM < R. This condition yields the maximum mass as 

M fermt = G- 3 / 4 m- 2 = m 3 pl /m 2 (1) 

where the Planck mass is defined as m p i = (ftc/G) 1 ^ 2 = 2.18 x 1CT 5 g. Thus the quantum 
mechanics (ft), gravity (G), and the particle physics (to) as well as relativity (c) characterize 
the universal structure of a fully degenerate fermion star (FDFS). 

Incidentally, degenerate bosons form much smaller structures since bosons have no exclu- 
sion principle like fermions. Therefore only the Heisenberg uncertainty principle (quantum 
pressure) can support the structure against the collapse due to gravity. The sole character- 
istic length scale is the Compton wave length Icompton = V ijnc) ~ R, which must be larger 
than the limiting scale, Schwarzschild radius, 2GM < R. This condition yields the density 
p = m 2 m 2 pl c 3 /ft 3 and 

M bose = m 2 pl /m, (2) 

which is known as the Kaup mass (Kaup 1968). This structure is a boson star. This structure 
is smaller than that formed by fermions by a factor m p i/m and no further discussion will be 
given in this paper, apart from a brief comment regarding the equation of state in §2. 

Now we proceed to consider the equilibrium structures made of fermionic dark matter. 
One may think that this problem is analogous to the white dwarf case and the Chandrasekhar 
mass is the limiting mass. However this is not the right answer. In the case of a white dwarf, 
the pressure is due to the degeneracy pressure of relativistic electrons, but the gravity is 
Newtonian because it is due to the rest mass of hadrons which are nonrelativistic. When 
we treat a degenerate star made purely of dark matter itself, the total rest mass is also 
due to dark matter particles themselves which can be relativistic. Therefore we need to 
consider the gravitational mass of the star in a general relativistic manner. This problem is 
therefore analogous to the case of a neutron star. In the case of the neutron star, there is a 
complication due to the fact that neutrons interact by nuclear forces. However in the case of 
weakly interacting fermionic dark matter, it can be considered as ideal gas and the equation 
of state is well defined. It turned out that this situation was treated by the classical paper 
by Oppenheimer & Volkoff (1939), since they assumed free neutrons and used the equation 
of state of ideal gas in their calculation. Basically fully degenerate fermions can form very 
compact dense structures including black holes (Bilic, Munyaneza & Viollier 1999; Bilic, 
Tupper & Viollier 2003). 

On the other hand, from the observational side, we now have many candidates of dense 
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structures at various scales. The most massive example is a huge dark matter distribution 
of the cluster of galaxies, A1689, reported recently by Broadhurst et al. (2005a,b). They 
obtained the mass column density distribution of the cluster of galaxies, A1689 by gravita- 
tional lensing. The column density profile has a flat-top and we suspected that this flat-top 
nature might be due to the degeneracy pressure of fermions. We derive a volume density 
profile from the observed column density profile assuming spherical symmetry and compare 
the observed 3D encircled mass profile with our model profile of an FDFS. 

The condense structure made from the fully degenerate fermion is not restricted to a 
center of a cluster of galaxies. If we consider more massive neutrinos, such as sterile neutrinos, 
the similar structures are realized in scaled-down form as in Eq.(l). If this structure is 
universal, we will find groups of black holes which have typical masses directly characterized 
by fermion masses. 

The paper is organized as follows. The formalism by Oppenheimer & Volkoff (1939) is 
introduced and equilibrium solutions are discussed in §2. Readers who are not interested 
in the derivation of the solutions, are advised to skip to §2.3 where the properties of the 
solutions are discussed. The application of a nonrelativistic FDFS to the cluster of galaxies 
A1689, is discussed in §3. The possible relation between the mass hierarchies of black holes 
and sterile neutrinos are considered in §4. 

2. Formalism 

Here we review the derivation of the general relativistic equilibrium equations following 
Tolman (1934) and Oppenheimer & Volkoff (1939). 

2.1. Relativistic treatment of equilibrium 

The most general static line element exhibiting spherical symmetry may be expressed 
in the form 

ds 2 = -e x dr 2 - r 2 d6 2 - r 2 sin 2 Odcj) 2 + e u dt 2 , X = A(r), v = u(r). (3) 

If the matter supports no traverse stresses and has no mass motion, then its energy 
momentum tenor is given by 
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T} = T* = T$ = -p, T 4 4 = e, (4) 

where p and e are respectively the pressure and the macroscopic energy density measured in 
proper coordinates. Einstein's field equations without the cosmological constant reduce to 
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(5) 



(6) 



dp 
dr 



p + e , 
— —v , 



(7) 



where primes denote differentiation with respect to r. These three equations together with 
the equation of state of the material e = p(e) determine the mechanical equilibrium of the 
matter distribution as well as the dependence of the metric g^s on r. 

The boundary of the matter distribution is the value of r = for which p = 0, and 
such that for r < r b ,p > 0. For r < r b the solution depends on the equation of state of the 
material connecting p and e. For many equations of state a sharp boundary exists with a 
finite value of ?v 

In the empty space, p = e = 0, surrounding the spherically symmetric distribution of 
matter, the Schwarzschild's exterior solution is obtained: 



e -A(r) = ! _ e u(r) = 1 _ ^ 



(8) 



where m is the Newtonian mass of the matter as calculated by a distant observer. 



Inside the boundary, Eqs. (5), (6), and (7) may be rewritten as follows. Using the 
equation of state e = e(p), Eq. (7) may be immediately integrated 



rp(r) 

v(r) = u(r b ) - / 



2dp 



p + e(p) ' 



e v(r) = e u(r b ) exp 



p(r) 



2dp 



o P + 



(9) 



The constant is determined by making e v continuous across the boundary. 



e w = (1 ) exp 

r 



D u(r) 



p(r) 



2dp 



o P + e (p) 



(10) 
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Thus e v is known as a function of r if p is known as a function of r. Further in Eq.(6). 
introduce a new variable 

1 2u 
u(r) = -r(l -e A ) or e A = 1 . (11) 

Then Eq. (5) becomes: 

d 4- = 4vre(p)r 2 . (12) 
dr- 
ill Eq. (5) we replace e _A by its expression (11) and v' by its expression (7). Then it 
becomes 

dp _ p + e(p) 3 



dr r(r — 2u) 



[47rpr s + u]. (13) 



Eqs. (12) and (13) form a system of two first-order equations in u and p. Starting 
with some initial values u — uq and p = po at r = 0, the two equations are integrated 
simultaneously to the value r = r& where p — 0, i.e., until the boundary of the matter 
distribution is reached. The value of u — Ub at r = r& determines the value of e A( - r ^ at the 
boundary to the exterior solution, making 



2m 



= m. (14) 



Thus the mass of this spherical distribution of matter as measured by a distant observer 
is given by the value Ub of u at r = r&. 

The following restrictions must be made on the choice of po and «o, the initial values of 
p and u at r = 0: 

(a) In accordance with its physical meaning as pressure, p > 0. 

(b) From Eq.(ll) it is seen that for all finite values of e~ A , u = 0. Since g n = — e A 
must never be positive, u < for infinite values of e A at the origin. However, it may be 
shown that of all the finite values of po at the origin p = is the only one compatible with 
a negative value of uq, and that for equations of state of the type occurring in this problem 
even this possibility is excluded, so that u must vanish. 

This can be seen from the following argument. Having chosen some particular value of 
p one may usually represent the equation of state in that pressure range by e = Cp s with 



-7- 



some appropriate value of s. Using this equation of state and taking the approximate form 
of Eq (13) near the origin for the case u < 0, and finite p , one obtains: 



Integration of this equation shows that for s < 1, p > can not be satisfied, and for 
s > 1 only the value po — is possible. This immediately excludes the possibility that 
degenerate bosons form an equilibrium structure, because p is independent of e, if p is solely 
due to thermal bosons (Landau & Lifshitz 1980). As we mentioned in §1, a boson star is 
supported by the quantum pressure of ground-state bosons, which is outside of the scope 
of the consideration of the equation of state (Kaup 1968). For the equations of state used 
for degenerate fermions, always s < 1 holds. It is also be noted that the above equation 
together with Eq. (10) show that e u ^ — > oo as r — > 0. 

(c) A special investigation for any particular equation of state must be made to see 
whether solutions exist in which < u < — oo and p — > oo as r — > 0. 



If the matter consists of fermions of rest mass /xo and statistical weight g, and their 
thermal energy and all forces between them are neglected, then a parametric form for the 
equation of state (Landau & Lifshitz 1980) is, 



dp p + e{p) p + Cp s 



(15) 



dr 2r 2r 



2.2. Equation of state for degenerate Fermi gas 



e = K(sinht — t), 



(16) 




(17) 



where 



K = 



8/i 3 



(18) 



and 



t = 4arcsinh , 

Hoc 



(19) 
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where pf is the maximum momentum in the Fermi distribution and is related to the proper 
particle density N/V by 

If we define the Fermi velocity t;? by 



Pf v f 



Hoc y/1 - iv F /c) 2 ' 



(21) 



t = 2 log ( 1 + VF ^ ) or Wc = tanh^. (22) 
\l-v F /cJ 4 

t and are independent of particle properties. 

Substituting the above expressions for p and e into Eqs. (12) and (13) one obtains: 

^ = 4vrr 2 ir(sinht-t), (23) 
dr 

and 



eft 



-4/[r(r - 2m)] 

x (sinh t - 2 sinh 1 /2t) / (cosh t - 4 cosh 1/2* + 3) 

4:71 

x — Kr 3 (sinht-8sinhl/2t + 3t) +m 



(24) 



These equations are to be integrated from the values u 
where t& = (which makes p — 0), and -u = 



0, t = t at r = to r 



So far, the equations are written in relativistic units, i.e., such that c = l,G = 1. This 
determines the unit of time and the unit of mass in terms of still arbitrary unit of length. 
The unit of length is now fixed by the requirement that K = 1/(4tt). From the dimensional 
analysis of Einstein's field equations, this requirement fixes the unit of length to be 



1 (2 

a = - - 
\g 



1/2 



3/2 



Hoc J {^Gy/ r 



(25) 
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and the unit of mass to be 




) 



1/2 



( 



h 



) 



3/2 



(M? 3 ) 1/2 ' 



(26) 



/i c 



Eqs. (12) and (24) written in a dimensionless form become: 



du 
dr 



r 2 (sinht — t), 



(27) 



dt 
dr 



-4/[r(r-2u)] 

x (sinh t - 2 sinh 1 /2t) (cosh t - 4 cosh 1 /2t + 3) 
x ir 3 (sinht-8sinhl/2t + 3t) +u . 



(28) 



For a given to, Eqs. (27) and (28) can be integrated. Fixing to is equivalent to fixing vpo- 
As long as dark matter particles are fully degenerate, the solution describes the equilibrium 
between the gravity and degeneracy pressure from Newtonian to general relativistic gravity 
and from nonrelativistic vf to relativistic vf- Therefore the solution for a given to (or ^fo) 
is independent of particle properties, while the units of length and mass (a and b) are fixed 
by the particle properties /x and g. 



The Eqs. (27) and (28) are numerically integrated using the fourth-order Runge-Kutta 
method. Apart from the gravitational mass u, there is another mass indicator, the dimen- 
sionless total rest mass y^, defined as 



which is the integral of the number density with the proper volume inside the radius tv 
to, vfo/c, lib, Ub, and r& are given for t = 0.1 ~ 14.0 in Table 1. At to = 0.1, degenerate 
particles are nonrelativistic {vfo/c = 0.025), while at to = 14, particles are extremely rela- 
tivistic (v F o = 0.998). In Fig. 1, the relation between the gravitational mass u b and outer 
radius r b is plotted. For t < 3, Uf, is an increasing function of t , while r b is a decreasing 
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Discussion on solutions 
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function of to. The maximum of u b = 0.0766 is reached for t = 3 and r = 0.663. This 
is the maximum stable solution for which vpo/c = 0.635, which is modestly relativistic. In 
Fig. 2, the radial profiles of the energy density, e, and pressure, p, are plotted. The contri- 
bution of p is relatively small compared to e. This solution is termed as quasi-Newtonian 
by Oppenheimer & Volkoff (1939). In the case of self-gravity of degenerate fermions them- 
selves, a stable configuration never reaches an extremely relativistic situation unlike the case 
of white dwarfs. Therefore the formula of the Chandrasekhar mass applied to this case 
over-estimates the maximum stable mass by one order of magnitude. For t > 3.0, there is 
no stable solution (Oppenheimer & Volkoff 1939). u b decreases takes a minimum value at 
(w fe ,r b ) = (0.0395,0.364) for t = 8.0. For a large t , (u b ,r b ) spirals to (~ 0.041, ~ 0.29). 

The gravitational mass defect, A = u b — y b is one measure of stability of a general 
relativistic equilibrium configuration. In Fig. 3, both y b and u b are plotted against t for 
0.1 < t < 14.0. For a small t , A is negative and gradually decreases and takes a minimum 
value, A = —0.0029 for the maximum mass stable solution (to = 3.0) and then increases. 
The fractional mass defect, or the packing fraction, is f — A/y b = —0.036. |/| is only 
3.6% and may appear small. However, if we compare this value to a typical nuclear packing 
fraction, which is less than 0.1% (Fermi 1950), we will find it to be significant. A < is a 
necessary condition for stability, but not sufficient. A < still holds for 3 < t < 5, but these 
solutions correspond to unstable equilibria. When A < 0, solutions are for stable equilibria 
for dA/dt < 0, while they are for unstable equilibria for dA/dt > 0. 

The question of whether dark matter fermions can form a black hole would be answered 
in terms of the maximum gravitational mass, M max and the corresponding radius R ma x, as 
functions of the particle mass /iq and statistical weight g. The results are: 



M, 



max 




(30) 



(31) 



and 



R. 



■max 




(kpc). 



(cm) 



(32) 



(33) 



It should be noted that Landau & Lifshitz (1980) define the maximum rest mass, M rest = 
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1.037 M max as the maximum mass, which is the total rest mass brought from infinity. In the 
case of a star made of free neutrons, fiQ = 939 MeV and g = 2. Then M max = O.7OM , 
and R max = 8.8 km. It should also be noted that R ma x is sensitive to the actual outer 
boundary condition in the numerical calculation, which in principle should be p = 0. In 
reality, integration must be stopped at a finite value of p and r b is somewhat sensitive to this 
finite p. Physically speaking, it is unrealistic to assume the fully degenerate equation of state 
to a small p and the outer boundary condition is not well defined. The total gravitational 
mass Ub is not very sensitive to the choice of the actual outer boundary condition p and 
is well defined. Although it is not explicitly stated in the original paper, Oppenheimer & 
Volkoff (1939) were probably aware of this limitation in the applicability of the equation of 
state, if we judge from the carefully chosen title "On Massive Neutron Cores" . We also admit 
that without this boundary condition, we cannot use the Schwarzschild's exterior solution 
outside the boundary and the formulation becomes more complicated. 

Before concluding this section, we comment on the nonrelativistic limit. For the non- 
relativistic limit, it is expected that between the mass, M, and the size, R, of an FDFS, the 
relation, 

MR 3 = const. (34) 

should hold (Landau & Lifshitz 1980). 

For t < 0.5, u b rf is indeed nearly constant and 

u b r 3 = 0.135, (35) 
within 4% precision. In physical units, 



MR 3 = (u b b)(r b a) 3 

= l-4x 10 124 ^ (^)" 8 (§-™ 3 ) (36) 

= 2.3xl0 26 l(^)" 8 (M -kpc 3 ). (37) 

Nonrelativistic solutions are similar and in Fig. 4, the profiles of the normalized mass 
density and 3D encircled mass are plotted. In the next section, we will find the behavior of 
the 3D encircled mass profile interesting and the logarithmic profile of the normalized 3D 
encircled mass is given in Fig. 5. 
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3. Can degenerate fermions be the dominant dark matter in the cluster of 

galaxies, A1689? 

Recently Broadhurst et al. (2005a,b) reported a mass column density profile of the clus- 
ter of galaxies, A1689, obtained from gravitational lensing. One of the important properties 
of the profile is that it has a flat top. We propose that this flat-top column density profile 
might be explained by the effects of degeneracy pressure of fermionic dark matter. Here we 
analyze this proposal. 

First we briefly introduce the main results of Broadhurst et al. (2005a,b). In their 
analysis, 1' corresponds to 129 kpc hr 1 . In Broadhurst et al. (2005a), the central 250 kpc 
hr x in radius of multi-color HST/ACS images were analyzed. The mass column density 
profile, S(r), is not expressed as a single power law of radius. The mass column density 
profile flattens toward the center with a mean slope of dlogS/dlogr ps —0.55 within r <250 
kpc h~ x . Inside the Einstein radius (9e ~ 50"), they obtained the slope of ~ —0.3 from the 
ratio between 9e and the radius of the radial critical curve, 9 r m 17". They fit their results 
with an inner region of an NFW profile (Navarro, Frenk & White 1996) with a relatively 
high concentration, C V i r = 8.2. 

The mass column density, S(r), is the integral of the volume density, p(r), along the 
line of sight over the entire cluster scale of Mpc. In order to study the possibility of fermion 
degeneracy near the center of the cluster, we need information on the volume density, p(r), 
instead of the column density, S(r). Broadhurst et al. (2005b) present the weak- lensing 
analysis of the wide field data obtained by Subaru and obtained the column density profile 
at r < 2 Mpc h' 1 . They fit the combined profile of HST/ACS and Subaru with an NFW 
profile with a very high concentration, C vir = 13.7, significantly larger than theoretically 
expected value of C„ r ~ 4. They also fit the same observed column density profile with a 
power law profile with a core. They give this result in terms of the angular radius dependence 
of the convergence, k, as 

Koc(e + e c )~ n . (38) 

9c = 1.65' and n = 3.16 give the best fit although 8c and n = 3.16 are mutually 
dependent and a finite range of the combination (9c, n) gives equally good fits. In terms of 
X 2 and the degrees of freedom, this core power law profile fits the observation better than 
the best-fit NFW profile and we use this profile for further discussion. Although Broadhurst 
et al. (2005b) do not claim so explicitly, the two facts that the best-fit NFW profile shows 
a much higher concentration than the value predicted by the CDM cosmology and the 
phenomenological profile, Eq.(38), fits better than the best- fit NFW profile, indicate some 
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serious contradiction to the CDM cosmology. 

We start our analysis from this core power-law profile, (38), for further discussion. We 
convert (38) to a column density profile, E(r), in physical units of length and mass using 
the relations, k = E/E CJ .j t , E cr jt ~ 0.95 g- cm -2 , and the normalization of 2D encircled mass 
inside the Einstein radius, t\e, f^ E T l (r)2nrdr = E^t • irr E . The result is expressed as 



E(r) = 25.2 • (r/r E + 2.2)~ 3 - 16 , (g • cm" 2 ) (39) 

where r E = 97 kpc /i" 1 corresponds to 9 E = 45", the value used in Broadhurst et 
al. (2005b). The 2D encircled mass, M2(r) = J E(r)27rr<ir, is analytically obtained and 
M 2 (r) = 1.3 x 1O 14 /i~ 2 M and 1.1 x 10 15 /i~ 2 M Q respectively for r = r E and r = oo. 
Therefore a high concentration of the mass is expected on the scale of r E . By assuming 
spherical symmetry, we wish to obtain the volume density p(r) by solving 



E(x) = J p(Vx 2 + z 2 )dz, (40) 

but we were not able to obtain an analytic solution. Instead, we assumed another power 
law with a core for p(r) and obtained the best fit parameters. The range of integration in z 
is from -2Mpc /i" 1 to +2Mpc h' 1 . The result is 



p(r) = 1.60 x 10" 23 (r/r E + 1.28)" 3 ' 71 h. (g • cm" 3 ) (41) 

Near the center, p(r) = 6.4 x 10~ 24 /i and 7.5 x 10" 25 /i (g- cm" 3 ) respectively at r = 
and r E . As for the case with the column density profile, the core radius and power-law 
index are mutually dependent and a finite range of their combination gives equally good 
fits. The above volume mass densities may appear small, but the degeneracy depends on the 
number density and de Broglie wavelength, both of which are heavily dependent on the rest 
mass of the particle. Before proceeding to the analysis by an FDFS, we first confirm that 
eV-mass fermions can be degenerate at these low mass densities. Since 1 eV corresponds 
to 1.8 x 10" 33 g, the number density, N/V xs 10 11 cm -3 , the mean inter-particle spacing, 
(N/V)~ l l z ~ 2 x 10" 4 cm. On the other hand, the de Broglie wavelength for a 1 eV particle 
with a velocity v is, h/pov = \com P ton-(c/v) = 2x 10" 5 (c/t> ) cm. Therefore for nonrelativistic 
particles with v < 0.1c, the condition for degeneracy, (iV/V)" 1 / 3 < A( dc Broglie), is satisfied. 

More quantitative discussion is possible for a 3D encircled mass profile given in Fig.6. 
For the purpose of an explicit comparison, here we fix h = 0.7. The 3D encircled mass profile 
also predicts the rotation curve profile as Fig. 7, which can be compared with observations of 
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kinematics of the galaxies in the cluster in the future. The column density profile, Eq.(39), 
was derived as a phenomenological formula without assuming any background physics, and 
so was the volume density profile, Eq.(41). Here we wish to explain the presence of the 
core in the power-law profile, which causes the flat-top nature of Eq.(39), in terms of an 
FDFS. The comparison must be made either in the volume density profiles or in the 3D 
encircled mass profiles. Here we choose the latter, because the observed 3D encircled mass 
profile is much less sensitive to the actual choice of the combination of the core radius and 
the power-law index in the original Eq.(39). As we mentioned above, the finite range of the 
parameter combinations gives equally good fits to the observations. The basic reason for 
this is that our observables are integrated quantities rather than the local densities. First, 
we compare the slopes of the 3D encircled mass profiles of the observation (Fig. 6) and our 
model (Fig. 5). The observed slope is 2.93 at r <10 kpc and 2.54 at 10 < r < 100 kpc. On 
the other hand, the slope of a nonrelativistic FDFS is 2.90 between r = 0.03r{, and 0.3r{, and 
becomes shallower at 0.3r& < r < r&. Note that the slope of the 3D encircled mass profile 
for a constant volume density is 3.0 and both the observed and model inner slopes of 2.9 
are consistent with the flat-top nature of the volume density profile. Now we proceed to the 
physical scaling of the model in terms of the mass and length. For this purpose, it is most 
convenient to use the relation Eq.(37) for the mass and size of a nonrelativistic FDFS. The 
Eq.(37) can be rewritten as 

logM 14 + 31o gj R 2 =f(g, /i ), (42) 

where M u = M/(1O 14 M ), R 2 = ^/(lOOkpc) and /(#, /i ) = 6.36-2 log g-S log(// /eV). 
In Fig. 8, the possible combinations of M and R are plotted for different values of f(g,no)- 
Intersections of the dotted lines and the solid line (3D encircled mass profile of A1689) corre- 
spond to solutions for representative values of /. Actually / is continuous and the solutions 
are continuous. The probable range is / = — 6 ~ 2. For this range of /, the rest mass 
of a fermion ranges from 30 to 2 eV. The larger the rest mass, the smaller the degenerate 
structure. The possible range of g is limited. For a particle with spin 1/2, g is either 1 
(Majorana particle) or 2 (Dirac particle). 

The special case is the massive neutrinos with similar masses for which effective g = 3 
(Majorana) or g = 6 (Dirac). Recent underground experiments have shown that the mass 
differences among three species of neutrinos are smaller than 0.05 eV (Shirai 2005). In order 
for massive neutrinos to have masses greater than 1 eV, they must have similar masses 
(degenerate in mass). The mean number density of relic neutrinos is cosmologically fixed 
and there is the well known relation for the contribution of the relic neutrinos to Q, 
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Q u h 2 = £^/(93.5eV), (43) 

where i — 1 ~ 3 for Majorana neutrinos and % — 1 ~ 6 for Dirac neutrinos. The 
maximum allowed neutrino mass can be estimated by setting Q u = 0.3 and h = 0.7, to be 
4.7 or 2.3 eV respectively for Majorana or Dirac neutrinos. The former gives / = or the 
degenerate mass of 1O 14 M in 100 kpc, and the latter gives / = 2 or that of 5 x 1O 14 M 
in 300 kpc. These are rather comfortable numbers for this cluster profile. If we question 
the CDM cosmology based on the possible inconsistency of the high concentration of the 
observed mass profile with the CDM predictions mentioned above, we might as well revive 
massive neutrinos (hot dark matter) as the candidate for the dominant dark matter. 



4. Mass hierarchy of black holes from that of neutrino through FDFS? 

In the previous section, we have studied a possibility that fully degenerate fermions forms 
a huge mass concentration at the center of a cluster of galaxies and we mainly examined the 
case for massive neutrinos. Is this possibility really true? In order to answer this question, 
we focus on the universality and the scalability of FDFS. We try to extend this idea of FDFS 
to other neutrinos and fermions with different masses. As we have already seen in Eq.(l,31), 
the rest mass of the fermion mostly determines the characteristic mass scale of the structure; 
more massive fermions form lighter structures. 

The most extreme condensed structure is a black hole. We have already known that 
there exist many black holes of several species in the Universe (Gebhardt et al. 2002). They 
are the most familiar stellar mass black holes (~ M Q ), giant black holes at the center of 
a galaxy (« 1O 7 M ), and the intermediate mass black holes (« 1O 3 M ). Although these 
massive black holes are actively studied based on the bottom-up scenarios that they are 
formed by the coalescence of the stellar sized black holes, we try to propose yet another 
scenario based on the context of FDFS. The most prominent property of those black holes 
are that they appear to have a hierarchy in mass range. Most of the black holes are classified 
in the above three types and those in other mass ranges is rare. Therefore it would be 
natural to suspect any definite mechanism to construct such hierarchy from the fundamental 
level. Our hypothesis is that such fundamental mechanism is the mass hierarchy in fermions 
or possibly neutrinos. If those black holes are formed by the overweight FDFSs, we can 
estimate the masses of those neutrinos as in Table 2. 

The corresponding masses of fermions, except v e ,n,T-, through Eq.(l) are far heavier than 
eV and we may identify those fermions as more massive sterile neutrinos. Lighter neutrino 
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of mass 10 _1 — 10~ 3 eV would yield FDFSs much extended dilute structures whose sizes well 
exceed the Horizon size. In the above, we have approximate values. However, the lower limit 
of the black hole mass within a class yields the precise value of the corresponding neutrino 
mass. The existence of such a lower limit would be the key ingredient of the FDFS model. 

5. Conclusions and Discussions 

We have examined the possible structures formed by the fully degenerate self-gravitating 
fermions (FDFS) at various scales, such as the mass concentration of a cluster of galaxies, 
the giant black holes at the center of a galaxy and the intermediate black holes in galaxies. 
As an order estimation, their characteristic masses directly reflect the constituent fermion 
masses through the simple relation, Mf erm i = G~ 3 ^m~ 2 = rrv^Jm 2 . For the purpose of 
a quantitative analysis, exact masses and detailed mass density profiles of FDFSs were 
examined from nonrelativistic to relativistic situation, using the formalism of Oppenheimer 
and Volkoff. 

These results were applied to the cluster of galaxies, A1689, whose mass distribution 
has been observationally obtained. We converted the observed column density profile to a 
volume density profile assuming spherical symmetry, and compared the observed and model 
3D encircled mass profiles. We found that the flat-top nature of the observed profile is 
reproduced by the model and the particle mass range is between 2 eV and 30 eV depending 
on the actual scale of the degenerate structure. 

For about black holes, our scenario will provide alternative mechanism of the black hole 
formation. Most of the present theories assume the coalescence of stellar mass black holes 
in the gravitational potential. Such processes seem to be quite complex compared to the 
FDFS scenario, and therefore it would be much difficult to explain, for example, the observed 
universal relation between the black hole mass at the center of a galaxy and the bulge mass: 
M BH /M bu i ss 0.002. This point will be discussed in detail in our future work. 

We thank Nobuo Arimoto for comments on the manuscript. This study was motivated 
by a very interesting talk given by Tom Broadhurst in 2004 at NAOJ, Mitaka, Japan. 
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0.02 0.04 0.06 0.08 0.1 

Gravitational Mass (u) 

Fig. 1. — Relation between the dimensionless gravitational mass u and the dimensionless 
radius r of equilibrium configuration for various central Fermi velocity (vf). For u < 0.04, 
there is only one r, while 0.04 < u < 0.0766, there are multiple values of r. Among them, 
however, only the configuration with the largest r is the stable equilibrium. The maximum 
dimensionless mass u = 0.0766 is the mass limit before collapsing to a black hole (square). 
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0.2 0.4 0.6 0.8 

Radius (r) 

Fig. 2. — Dimensionless energy density (e) and pressure (p) as functions of dimensionless 
radius for t = 3.0. This is the maximum mass stable configuration, p is small compared to 
e. Degenerate particles are moderately relativistic and the gravity is quasi-Newtonian. 
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Fig. 3. — Dimensionless gravitational mass (ub) and total rest mass plotted as functions 
of to- The gravitational mass defect, A = Ub — Ub, is negative for t Q < 3 and takes the 
minimum value for the maximum mass stable configuration (to = 3). 
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0.2 0.4 0.6 0.8 

r / r b 

Fig. 4. — Normalized density and encircled mass profiles for a nonrelativistic FDFS. Non- 
relativistic solutions are similar. 
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log(r/r b ) 



Fig. 5. — Normalized 3D encircled mass profile of a nonrelativistic FDFS in logarithmic 
scale. The slope of the profile is 2.9 for log(r/r b ) = —1.5 ~ —0.5. 
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Fig. 6. — 3D encircled mass profile of Af689 obtained from the best-fit volume density 
profile. h=0.7 is assumed. 
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Fig. 7. — Predicted rotation curve of A1689 estimated from the 3D encircled mass profile. 
h=0.7 is assumed. 
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log(r) (kpc) 

Fig. 8. — Allowed combinations of the mass and size of an FDFS for different values of 
f(g,Ho) = 6.36 — 21ogg — 81og(/i /eV) (dotted lines). The solid line represents the 3D 
encircled mass profile of A1689. Intersections between the solid and dotted lines correspond 
to possible degenerate structures. 
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Table 1: Solutions for various t . 



to 


vfo/c 


u b 


Vb 


r 




Fermi Vel. 


Grav. mass 


Rest Mass 


Radius 


0.1 


0.025 


0.0012 


0.0012 


4.934 


0.2 


0.050 


0.0033 


0.0033 


3.484 


0.3 


0.075 


0.0060 


0.0060 


2.838 


0.4 


0.100 


0.0091 


0.0091 


2.450 


0.5 


0.124 


0.0126 


0.0126 


2.183 


1.0 


0.245 


0.0325 


0.0328 


1.495 


1.5 


0.358 


0.0516 


0.0525 


1.161 


2.0 


0.462 


0.0659 


0.0677 


0.943 


2.5 


0.555 


0.0740 


0.0766 


0.784 


3.0 a 


0.635 


0.0766 


0.0795 


0.663 


4.0 


0.762 


0.0710 


0.0730 


0.493 


5.0 


0.848 


0.0598 


0.0597 


0.391 


6.0 


0.905 


0.0491 


0.0470 


0.334 


7.0 


0.941 


0.0420 


0.0387 


0.343 


8.0 


0.964 


0.0395 


0.0360 


0.364 


9.0 


0.978 


0.0416 


0.0383 


0.382 


10.0 


0.987 


0.0444 


0.0414 


0.367 


11.0 


0.992 


0.0471 


0.0443 


0.399 


12.0 


0.995 


0.0463 


0.0435 


0.362 


13.0 


0.997 


0.0442 


0.0411 


0.324 


14.0 


0.998 


0.0411 


0.0376 


0.285 



"This solution corresponds to the maximum mass stable configuration. The gravitational mass defect, Ub — yb, 
takes the minimum. 



-27- 



Table 2: Neutrino masses and possible hierarchy 





m u /eV < 


Mf erm i/M Q 




FDFS 




2.3 


9.8 x 10 16 


42 kpc 


the center of a cluster 


V\ 


0.18 x 10 6 


2.5 x 10 7 


460^ 


giant BH at center of galaxy 


V2 


18.2 x 10 6 


2.7 x 10 3 


0.050-R© 


intermediate BH in a galaxy 



